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We study an “observable-based” notion of equilibration and its application to realistic systems 
like spin qubits in quantum dots. On the basis of the so-called distinguishability, we analytically 
derive general equilibration bounds, which we relate to the standard deviation of the fluctuations 
of the corresponding observable. Subsequently, we apply these ideas to the central spin model 
describing the spin physics in quantum dots. We probe our bounds by analyzing the spin dynamics 
induced by the hyperfine interaction between the electron spin and the nuclear spins using exact 
diagonalization. Interestingly, even small numbers of nuclear spins as found in carbon or silicon 
based quantum dots are sufficient to significantly equilibrate the electron spin. 

PACS numbers: 03.65.Yz, 05.30.-d, 76.20.+q, 85.35.-p 


I. INTRODUCTION 

The theoretical understanding of the notion of equili¬ 
bration in closed quantum systems has significantly de¬ 
veloped in recent years.®^ In the absence of a thermal 
bath and the presence of quantum fluctuations, the classi¬ 
cal concepts of the physical and mathematical description 
of equilibration do in general not work anymore. 10 ^ There¬ 
fore, it is of utmost importance to first properly define 
what we mean by equilibration in closed quantum sys¬ 
tems that can even be in a highly non-thermal state. This 
difficult task is one of the driving forces of the research 
area of quantum thermodynamics. Useful concepts im¬ 
ply different definitions of equilibration. For instance, 
many authors identify equilibrium with the saturation 
of the expectation values of certain observables 1SU2H 12 ] 
These ideas are appealing as they are intuitive and the 
relevant quantities are measureable. However, as it is ar¬ 
gued in Ref. El this definition is not satisfying because the 
measurable probabilities of the outcomes of an observ¬ 
able may still be dynamical while its expectation values 
has saturated. In our article, we start from a more so¬ 
phisticated concept^ and link it to the above discussed 
ideas of saturating expectation values. Doing so, we con¬ 
nect abstract definitions of quantum equilibration to a 
very concrete experimental system where one could po¬ 
tentially see this exciting physics. 

The system we have in mind is an electron spin 
confined in quantum dot (QD) that functions as a 
spin qubitIn typical host materials, like GaAspHUl 
silicorP-dL22J or carbonpSHSD e i ec ^ ron S pi n j s co upled 
to many nuclear spins by the hyperfine interaction. 40 41 
The number of nuclear spins that matters for the elec¬ 
tron spin depends on the size of the QD, i.e. the confine¬ 
ment, and the natural abundance of nuclear-spin carrying 
isotopes of the host material. In practice, this number 
can range between very few nuclear spins (in silicon- or 
carbon-based systems (paMI to millions (in GaAs-based 
QDs)!®^ In recent years, the experimental control of 
spin qubits in QDs has developed to a state of perfection 
at the single- and two-qubit levelP^D^It is possible to ini¬ 


tialize, manipulate l 22 l 49 t ESl and read-oulPSESI S pj n qubits 
with a very high precision and to even engineer the state 
of the nuclear spin battP^-H^to increase coherence times. 
These remarkable experimental achievements have been 
accompanied by sophisticated theoretical research®"^ 
which has promoted the understanding of spin dynamics 
in QDs with regard to the hyperfine interaction. There¬ 
fore, we believe that these systems are ideally suited to 
study predictions related to quantum equilibration. 

First, we have to introduce a general theory of equi¬ 
libration of the closed quantum system that fulfills all 
the requirements of the realization that we have in mind. 
This will be done on the basis of the distinguishability ^ 
which is a measure to distinguish the actual state of 
a quantum system from its equilibrium state on the 
basis of a finite set of observables. If the values of 
the distinguishability are on average smaller than a 
given reference value e we argue that the quantum sys¬ 
tem is e-equilibrated. In order to connect our concept 
of equilibration with experimentally measurable predic¬ 
tions, we first rel ate t he distinguishability with the weak 
distinguishability p 5 8 which offers an equivalent descrip¬ 
tion of equilibration under special conditions (for two- 
outcome observables). The time-averaged weak distin¬ 
guishability (TAWD), however, is capable to bound vari¬ 
ances of expectation values from above. We have an¬ 
alytically derived certain bounds for the TAWD, which 
depend on the Hamiltonian and the initial state of the 
quantum system. As a consequence, our analytical equi¬ 
libration bounds for the TAWD should directly affect the 
experimentally determined variance of the measurement 
operator. Therefore, it should be possible to modify the 
system at hand such that the bounds are varied and to 
see the difference in a direct measurement of the variance. 
Evidently, this is a concrete prediction of an observable 
consequence of quantum equilibration. 

With this prediction at hand, we eventually try to bet¬ 
ter understand quantum equilibration by looking at our 
central spin model mentioned above. In order to cal¬ 
culate the TAWD here, we treat very simple observables 
like the electron spin operator in direction parallel or per- 
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pendicular to an external magne tic field. Since we em¬ 
ploy exact diagonalizatiorpSEaT 87 dSll for this calculation, 
we are limited to a finite number of nuclear spins (up to 
10). However, in state of the art QDs based on silicon 
or carbon host materials, such numbers of nuclear spins 
are within experimental reachP^HH Hence, the finite size 
effects we analyze in this article should eventually be ex¬ 
perimentally relevant. We find that our analytical re¬ 
sults of bounds of the quantum equilibration describe 
very well the numerical simulations based on the central 
spin model for compatible conditions. This makes us 
confident that our predictions can be really seen in mea¬ 
surements of the spin dynamics of a confined electron 
spin coupled to a bath of nuclear spins. 

The article is organized as follows. In Sec. II, we ex¬ 
plain the notion of equilibration employed in this work 
and introduce the (weak) distinguishability used to de¬ 
scribe it. Subsequently, in Sec. Ill, we will derive ana¬ 
lytical results of equilibration bounds. In Sec. IV, these 
general results are then compared to a central spin model 
of an electron spin in a QD coupled to a quantum bath 
of nuclear spins. We conclude in Sec. V with a summary 
of our main results. Some derivations are presented in 
three Appendices. 

II. BASIC CONCEPTS OF EQUILIBRATION 

In this section, we briefly describe known concepts of 
quantum equilibration for future reference. We consider 
a closed quantum system whose state pit) evolves ac¬ 
cording to the von Neumann equation p(t) = 7 [p(t),H] 
where H is the d dimensional Hamiltonian of the total 
system TL. Due to the unitary time evolution, each fi¬ 
nite quantum system obeys a recurrence time Tr > 0, at 
which the state of the system approaches within some ac¬ 
curacy its initial state. However, this time does not play 
a role in most experiments as it scales exponentially— 
with the dimension of T~L and is almost always much larger 
than the age of the universe. With the commonly used 
and well-defined time-averaged statd^ 

w := (p )oo = lim (p) t (1) 

£—>■00 

one circumvents recurrence problems. Throughout this 
article, (/) t = \ dt' fit') is used to denote time aver¬ 
ages. This time-averaged state can be considered as an 
equilibrium state for several reasons. First, it does not 
evolve in time as [cD, H] = 0. More importantly, if the 
expectation value O(f) := tr [Op{t)\ of any observable O 
saturates at some value for long times z it can be calcu¬ 
lated by (O )oo = tr[(p(t)) oc O] = tr [CjO\. In contrast to 
thermal states like the Gibbs state, Co generally depends 
on the initial state p. 

Analogously to earlier work^HSl we regard a quantum 
mechanical system to be in equilibrium if one cannot dis¬ 
tinguish between the state p(t) of the full system and its 
equilibrium state Cj for most times by applying a finite 


set of measurements T = {Oi} that can be performed in 
an experiment. These measurements are not restricted to 
subspaces of the whole Hilbert space. Hence, this defini¬ 
tion does not rely on the subdivision of the full quantum 
system into a small, measurable system and a large, not 
measurable bath. 

For the above notion of equilibration, it is not sufficient 
that the expectation value of an observable O = JA XjPj 
saturates, since p(t) and lo can still be distinguished by 
the (experimentally) measurable probabilities tr [p(t)Pj] 
of its eigenvalue A j. Rather each of these time-depend 
probabilities has to saturate in order to guarantee in- 
distinguishability. Considering this necessity, SliortP has 
introduced the distinguishability 

djripit), Co) = max ^ ^ | tr^-p)*)] - tr [PjCj] \. (2) 

as a proper measure of distance between pit) and Cj. 
Mathematically, it is closely related to the trace distance, 
but considers the finite number of accessible measure¬ 
ment operators. In contrast to the trace distance, how¬ 
ever, this measure is not a metric, but a semi-metric since 
djripit), Cj) = 0 is possible for p(t) jr Co. This behavior 
is important, because it permits the desired property of 
equilibrated states: A sufficient condition for equilibrium 
is that one is not capable to distinguish the state of the 
system pit) from Co for most times by the set CF of mea¬ 
surements. 

In order to account for the fact that the state of the sys¬ 
tem must be indistinguishable for most times during the 
time evolution, one can demand the time-average of the 
positive quantity djripit), Cj) to be smallPConsequently, 
we regard a system to be e-equilibrated at time t if 

(' d^ipit),u))) t < e , (3) 

where £ is a small positive constant, which we are free to 
choose. A reasonable choice for this constant is, for in¬ 
stance, the precision of the measurement devices in an 
experiment. Further, we call systems equilibrating in 
a time interval / if the time-averaged distinguishability 
( djrip,Ci))t decreases on average within I. 

By introducing the distinguishability and its time av¬ 
erage in Eq. (|3|, we achieved a suitable mathemati¬ 
cal definition of our concept of equilibration. How¬ 
ever, the distinguishability cannot be measured directly 
in an experiment. Yet, with a slight modification of 
the distinguishabil ity, one can find the so-called weak 

distinguishabilitjl^MlfilSl 

D d it) := (tr [p{t)0] - tr[wO]) . (4) 

This quantity is unlike the distinguishability djripit),Co) 
only given by the expectation values of C® with respect 
to both p and Cj, but does not depend on the probabili¬ 
ties to measure individual eigenvalues. Hence, the weak 
distinguishability carries the same unit as the squared 
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FIG. 1: (Color online) Connection between equilibrium, dis- 
tinguishability, weak distinguishability, and expectation val¬ 
ues. For clarity within this figure, we assume that only one 
observable O is measurable, i. e. T = {O}. If O is a measure¬ 
ment with only two possible outcomes a ^ b, the weak distin¬ 
guishability is equivalent to the distinguishability. Therefore, 
saturated expectation values are then a necessary and suffi¬ 
cient condition for equilibrium. 


measurement operator and takes values between 0 and 
4||0|| 2 , where ||0|| is the spectral nornP^ of O. As we 
will show below, the long-time average of the weak dis¬ 
tinguishability can be identified with the variances of 
the observable, which can be determined in an exper¬ 
iment. A small time-averaged weak distinguishability 
(TAWD) (D§) t < e' <C ||0|| 2 is a necessary condition for 
(do(p(t),oj)) t < £ <C 1 and, hence, according to Eq. ^ 
for the system to be in equilibrium. If two-outcome mea¬ 
surements O = aP a + bPb (a ^ b,P a + Pb = Id) are 
considered, both quantities are even equivalent as they 
are then related to each other by 

D d (t) = ( a ~b) 2 [d d (p(t),u>)] 2 . (5) 

In Fig. [T] we summarize these dependencies and the con¬ 
nection to equilibrium. As a last property of the TAWD, 
we show in App. [A] that the TAWD 


( 6 ) 

is separable in a time-independent part A^ and a time- 
dependent part Sq(J:), which decreases at least with 
8&{£) = O (f _1 ) for t —> oo. This behavior will play 
an important role for relating the TAWD to measurable 
quantities in the next section. 


III. EQUILIBRATION BOUNDS 
Weak distinguishability vs. variance 

As argued above, the TAWD is a useful quantity to 
describe equilibration in closed quantum systems. More¬ 
over, it is directly related to measurable properties of the 
system under consideration. As we explicitly derive in 
App. E the variance Var^f, At) of expectation values 


0{t') in a time interval t' G I = [t, t + At] is bounded by 

Var 6 (t,Af)< (D d ) t , (7) 

where the size At of the time interval I needs to be suf¬ 
ficiently large. More precisely, At must be of such a size, 
that (Do) t does not increase on average within /. The 
above estimate even turns into an equality if (D^t is con¬ 
stant within I. According to Eq. © , this is the case for 
each system and all observables at long times because the 
TAWD converges. Consequently, its infinite-time limit 

Aa = lim {Dp.) t = lim VarMt, At) (8) 

equals the variance of expectation values in any time in¬ 
terval / at long times. Since the time-dependent part 
S^(t) can decay much faster than 0(f -1 ), this satura¬ 
tion will be already reached within finite times for many 
systems. Note that A^ is the variance of expectation 
values of O that does not arise from measurement errors 
but from fluctuations within the finite quantum system. 
Hence, a valid interpretation of A q is a measure of the 
capability of the system to equilibrate with respect to O. 
The smaller A^ the less fluctuations of the expectation 
values of 0(t ) around tr[6jO\ are present. 


Useful equilibration bounds at large times 

As we elaborately show in App. [C] the long-time val¬ 
ues of the TAWD can be estimated in different manners 
giving rise to the bounds 

Aq < A^ := 7Vg||/ 5|| 2 tr[0 2 ], (9) 

Aq < A^ := iV G ||0|| 2 tr[/5 2 ], (10) 

A 6 <Ay=« G Hl. (id 

Before we discuss and compare these findings, we focus 
on the quantities they depend on. First, in all bounds the 
maximum degeneracy Nq of gaps in the energy spectrum 
of the Hamiltonian enters, whose size is, hence, crucial for 
them to be of reasonable magnitude. Note that it is not 
sufficient to have a non- deg enerate eigenvalue spectrum 
in order to reach Nq = 1IL22I The properties of the observ¬ 
able enter the equations by ||0|| and tr[0 2 ], which are re¬ 
lated to each other by ||0|| 2 < tr[0 2 ] < ||0|| 2 rank 6 . 101 
The bounds also respect the consequences of different 
initial states p. Explicitly, tr[p 2 ] is its purity and ||/5|| 
is the maximum eigenvalue of the initial state, where 
tr[/5 2 ], ||p|| > d -1 . Moreover, the initial state also de¬ 
termines the size of the so-called effective dimensiorP 
d e g, which is defined by d~g = JT(tr[.Ej/3]) 2 with Ej 
being the projector onto the eigenspace of energy E r 
Vividly, d e g quantifies the dimension of the Hilbert space 
that is actually reached during the time evolution. It 
reaches values between 1 and d. The latter is the case 
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for the totally mixed state p = ^ 1 or for pure states like 

Due to the last property of d e g, the third estimate 
Ajl, which has previously been found in Ref. (H with a 
different approach, is the most restrictive bound if pure 
initial states are considered. The other two estimates are 
useful for mixed states as both ||/5|| and tr[p 2 ] become 
small if and only if the state p is mixed. The advan¬ 
tage of is that the quantity tr[p 2 ] is independent of 
the basis whereas one needs to know all eigenstates and 
eigenvalues of H in order to calculate d e g. The bound Ab 
is more restrictive than Ab, Ab, and previously found 
estimate:!^, if 


d eff ||/3|| 2 rankd < 1. (12) 


Thus, the rank of O should be small while the mixture 
of the initial state p should be high, since d e g scales as d 
for very mixed states while ||p|| 2 scales as d~ 2 . 


Generalization to finite times 


So far we have focused on the behavior of the TAWD 
for long times. However, according to Eqs. (|6j) and Q, 
we can even give estimates for finite times provided that 
one can bound the time-dependent part 5q (t). As we 
have discussed, 5@(t) is at least decaying as f” 1 in the 
long-time limit. In a recent analysis, L. P. Garcfa-Pintos 
and coworker^Sl have derived many interesting proper¬ 
ties of the TAWD. Among other things, the authors 
have bounded the time dependent part of the TAWD 
by 8q (t) < j. where L is a constant that dependents on 
p,H, and O. Combining this equation with our previ¬ 
ous results, we find: Given a system with arbitrary ini¬ 
tial state /?, Hamiltonian H and an arbitrary observable 
O, we can estimate the variance of expectation values 
around the long-time average tr[wO] within any time in¬ 
terval [ t , t + At] by 


Vai d (t,At) < ( D d ) t < 


L 

1 


+ min A*x. 
i ^ 


(13) 


The infinite-time bounds A*^ are given in Eqs. (|9| to ( jll[ ) . 

If one of the A*^ turns out to be a small number and O is 
a two-outcome measurement, the system will equilibrate 
in the sense defined in Sec. II In that sense, A*^ gives 
an estimate for the ability of a closed quantum system 
to equilibrate. Even if this concept of equilibration is 
not used, the above bounds still estimate the variances 
of observables in any closed system correctly. 


IV. APPLICATION TO SPIN MODELS 


A. Central spin model basics 


In this section, we apply the general concepts of equi¬ 
libration explained above to a specific, realistic system. 
This allows us to show the physical significance of the 
above ideas for experiments. In particular, we make con¬ 
crete prediction s on measurable properties of an electron 
spin in a QD ,E!3 I l6 l40|4i | w ] 1 j c } 1 j s interacting with the nu¬ 
clear spins of the host material. 

Besides this hyperhne interaction (HI) between the 
electron spin and the nuclear spins, we consider an ex¬ 
ternal magnetic field, which is commonly used to split 
the Zeeman levels of the spins. In many experimen¬ 
tal setups, further effects such as direct interactions be¬ 
tween nucle ar sp ins and spin-orbit mediated effects are 
negligible]- 14 15:8 -' These interactions are, thus, not taken 
into account in our model. By this choice of the inter¬ 
actions, we construct a minimal model, which is realized 
by several experimental setups jm t 16 l 40 l 41 l 

Prominent examples are devices made from group IV 
elements, which exhibit nuclear spin-less isotopes. Iso¬ 
topic purification of carbon or silicon allows to manipu¬ 
late the number of nuclear spins present in the QDP^IEU 
This possibility allows to probe the influence of the sys¬ 
tem size on our bounds in Eqs. ^ to ((nj. 

In the following, we are especially interested in how the 
nuclear spins will equilibrate the electron spin. Since the 
observables of the electrons spin S XtViZ all have two out¬ 
comes, the distinguishability and the weak distinguisha- 
bility are equivalent according to Eq. ([5]) . The saturation 
of the expectation values of spin operators, hence, corre¬ 
sponds to the equilibration of the full system - given that 
they are the only accessible measurements. 

After these general considerations, let us introduce the 
total Hamiltonian H = Hjjj + Hze describing our model 
in more detail. Although our qualitative results areinde- 
pendent of this choice, we choose a graphene Ql j 89 l 91 l lMJ 
as a reference in order to benefit from previous results 
Then, the HI Hamiltonian is given by 


Hm = Am J2 \fa\ 2 M - \^+i- + S-i k + )\, (14) 

k =1 

where the energy scale of the Hi-—is Ahi = 0.6 peV 
and the number of nuclear spins is K. We use dimen¬ 
sionless spin operators S XtVjZ , I x ,y,z, S± = S x ± i S y , and 
I±, analogously. The probability to find the electron at 
the site of the fc-th nuclear spin is given by the abso¬ 
lute value of the envelope function |<^fc| 2 > \(f>k-i\ 2 - The 
strongest HI coupling Aml^ifl 2 defines the characteristic 
time thi = h/(Am\cj ) K | 2 ). Whenever we average over dif¬ 
ferent initial conditions, we maintain a maximum ratio 
of |</>fc| 2 /|</>j| 2 < 100 for all k,j. For qualitative results, 
we present the results for an exemplary set of coupling 
constants as we have found similar results for many ran¬ 
domly generated sets of coupling constants. 













5 


The effect of an external magnetic field B z is described 
by the Zeeman Hamiltonian 

Hze = b ■ A m \<f> K \ 2 f S z - 7 4* j , (15) 

where b is the strength of the field in units of the strongest 

HI coupling. Note that the nuclear spins couple only very 

weak to external magnetic fields compared to the electron 

spin: in carbon we find 7 = gNg,JV w4x 10 -4 . 

. . gm 

Besides the Hamiltonian, the time evolution of ob¬ 
servables depends on the initial state p. In the follow¬ 
ing, we choose product states p = p e \ ® p nuc , where 
Pei = \'4 , ei){4’ei\ and p nuc describe the uncorelletad ini¬ 
tial states of the electron and nuclear spins, respectively. 
This assumption is plausible since the initial state of the 
electron spin can be experimentally well prepared in a 
pure (polarized) state by means of an ext ernal magnetic 
fielcfST using light in optically active CjD d 1 5 I 5 S ]ES 1 p r Py 
suitable pulse sequences in double QD setups LL 6 l 22 l 49 t t57] 
The nuclear spins, however, will on average be in an un¬ 
polarized state if no further efforts are undertaken in an 
experiment. Since experiment ally relevant temperatures 
are on the order of ml( to K! 14 l 15 l, the thermal energy 
exceeds all other energy scales of the nuclear spins by 
far. On top of that, to follow the time evolution of the 
electron spin, many repetitions of the experiment are 
needed. Since each of these runs start with a different 
initial state, the nuclear spin state can be described by 
a totally mixed state p nuc = 1/2 K on average. However, 
the nuclear spin state can be also mani pul ated by means 

of dynamical nuclear polarizatioii^SMMIZO <i anc ] state 

narrowing ) 64 * 69 * 7 2-75 which allow to significantly polarize 
the nuclear spins and to change the composition of the 
initial state of the nuclear spinsPA^ Motivated by these 
experimental possibilities, we also investigate the effect 
of polarized initial states of the nuclear spins by using a 
Gaussian distribution of states. 

With both, the Hamiltonian and the initial state given, 
the time evolution of the density matrix and, hence, of 
every observabl e in the system can be calculated by exact 

diagonalizatioiPSlSS sTHU], which j s performed using the 

EIGElsP package for C++. 


B. Spin dynamics 

Once the time evolution of an observable is known, 
its variance, the weak distinguishability and the TAWD 
defined in Sec. |TT] are readily calculated. This enables 
us to demonstrate, that the TAWD indeed bounds the 
variances of an observable. 

As an example, we show the evolution of S x (t) in 
Fig- i At times t ~ thi, the initially polarized elec¬ 
tron spin begins to oscillate with decreasing amplitude 
around its long-time average (S ' a; ) 00 ~ 0. The square 
root of the weak distinguishability \J(Dg ) t bounds the 
standard deviation of S x (f) as predicted at all times. At 



FIG. 2: (Color online) The time evolution of the electron spin 
component S x (t) and its long-time average tr[cbSy. We use 
a magnetic field b = \ and K = 6 nuclear spins with random 
coupling constants. Initially, the electron spin is maximally 
polarized in ^-direction while the nuclear spins are in the to- 
tally mixed state. The square root of the TAWD V(D§ )t 
bounds the standard deviation of S x (t) for all times and con¬ 
verges to it for large times. For comparison, we plotted the 
analytically derived bounds given in Eq. (131 and Eq. (101. 
Note that all quantities are dimensionless. 


large times, the TAWD (Dg } t saturates to a finite value 
whose size corresponds to the quantum fluctuations in 
our finite model. As explained above, the TAWD in turn 
can be bounded itself by the analytical expression given 
in Eq. (13). For finite times, this bound deca ys w ith 
0(t _1 ), while it saturates at A| given in Eq. ( jioj ) for 
large times. For the parameters chosen in Fig. [2j we find 
A| = Remarkably, already for K = 6 nu- 

clear spins, this long-time estimate yields a very sharp 
upper bound on the standard deviation of fluctuations of 
the signal. 


As explained above, the properties of the TAWD and 
its bounds depend on the Hamiltonian of the system. 
Thus, one should test how different Hamiltonians alter 
the equilibration. In a QD, the easiest way to change the 
Hamiltonian is to modify the external magnetic field. By 
varying b over approximately two orders of magnitude, we 
sweep from a situation in which the electron spin couples 
most strongly to the nuclear spins to a scenario where the 
Zeeman coupling is dominant. In Fig. |3j we compare the 
TAWDs of S x and S z . For both spin components, we ob¬ 
serve that equilibration sets in approximately at time thi 
and reduces the initial values of the TAWDs roughly by 
two orders of magnitude for all values of b. As we discuss 
later, the size of this reduction depends on the number 
of nuclear spins. In fact, even high values of b cause only 
Lamor oscillations of S x at small times tze oc b, but do 
not change the overall equilibration behavior. Besides 
this, the only effect of large magnetic fields is a reduced 
initial value of the TAWD for S z . This can be under- 
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t [thi] 

FIG. 3: (Color online) The TAWD as a function of time for 
different external magnetic fields b for a) S x and b) S z . The 
electron spin is initially polarized in x- and 2 -direction, re¬ 
spectively. Apart from that, the same parameters as in Fig. [2] 
have been used, a) All TAWD s start with the same value 
and show oscillations on a timescale tze oc b. The actual 
equilibration starts at thi- b) The TAWDs start with a dif¬ 
ferent value for different b. The equilibration due to the HI 
also starts at thi. 


stood as follows. As the electron spin is initially fully 
polarized parallel to a strong magnetic field, its initial 
state is almost fully preserved, since the flip-flop terms 
of the HI are suppressed due to the large Zeeman splitting 
of the electron spin states. In other words, the electron 
is initially approximately in an eigenstate of the total 
Hamiltonian for strong external magnetic fields. Hence, 
p is already initially close to ui and, as a consequence, 
indistinguishable from Gj by means of S z . 

Similar effects can be found for polarized states, which 
we approximate by a Gaussian distribution of states char¬ 
acterized by a mean polarization p and a standard de¬ 
viation Cp. First, we have fixed a p = 0.3 and varied 
the mean polarization between p = 0 and p = 0.75 for 
a system containing K = 6 nuclear spins. Since the 
initial states approach eigenstates of the total Hamilto¬ 
nian for increasing polarization, we observe a decreas¬ 
ing initial distinguishability Dg (0). Due to the HI, the 
TAWD (Dg ) OQ saturates again around a value, which 
is about two orders of magnitude smaller than its ini¬ 
tial size. For larger polarizations this reduction becomes 
smaller, since the HI spin flip-flops become less effective. 
These findings are consistent with a smaller effective di¬ 
mension d e ff of polarized initial states in Eq. ©• Anal¬ 
ogous simulations with standard deviations in an interval 
0.15 < a p < 0.75 show no significant differences to these 
observations. 


C. Size dependence of the bath 


We finally want to address the question how many 
nuclear spins are required in order to treat them as a 
bath. By adding more and more nuclear spins, no sud¬ 
den change is observed but the fluctuation of spin com¬ 
ponents of the electron decrease exponentially with the 
number of nuclear spins, cf. Fig. [4] 

The numerically obtained values of the long-time 
TAWD are about one order of magnitude smaller than 
the presented bound A| . Considering the estimates 
made and the generality of the bound, this is still a fairly 
good result. Fig. [3] also suggest that quantum fluctu¬ 
ations may decrease even faster with increasing system 
size than our analytic bounds require. Note that this 
A'-dependence of the equilibration properties is not lim¬ 
ited to mixed states only. Reconsidering previously ob¬ 
tained datrPil, we have calculated the effective dimen¬ 
sion for randomly chosen pure initial states of the nu¬ 
clear spins. For these states, it scales approximately with 
d e ff ~ d/2 = 2 K . According to A| oc d~g = 2~ K given 
in Eq. ( [IT] ) , this dependence also gives rise to an ex¬ 
ponential decay of (Dg ) oc , which is confirmed by our 
numerics As discussed by ReimamP, the effective di¬ 
mension of almost all states grow exponentially with the 
size of the system. Hence, such a decay is a rather generic 
result, which can be understood as follows. If we add a 
nuclear spin to the system, we double both the size of 
the Hilbert space and the number of energies driving the 
dynamics of the electron spin, which finally leads to the 
observed reduction of fluctuations. 

We can indeed generalize these findings to other quan¬ 
tum systems that differ from our model, e.g. a central 
spin model with isotropic hyperfine interaction or even 
topologically different models like spin chains. Given 
that the effective dimension d e e ~ d scales exponentially 
with the bath size, either due to totally mixed bath states 
or due to randomly chosen pure initial states,^ we can use 
the bound in Eq. © to deduce the following statement. 
The number N of bath spins that are sufficient to satu¬ 
rate an electron spin in some arbitrary quantum model 
increases only logarithmic with the inverse resolution 1/r 
of the measurement: 


N > log 2 


^dlOII 2 



(16) 


where c = d e s/d < 1. For a resolution r = 0.01 H and an 
initial state far away from an energy eigenstate (c —I 1), 
the electron spin components equilibrate in any quantum 
model with non-degenerate gaps (Ng = 1) if the electron 
is coupled to more than 11 bath spins. As our model 
demonstrates, even less bath spins iV « 7 are capable 
of equilibrating the electron spin components below this 
resolution in experimentally relevant scenarios. 
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FIG. 4: Color online) Dependence of the long-time TAWD 
for the observable S z on the number of nuclear spins K. We 
approximate (D §_)oo by (D§ )tei, where I = [r,r + At] is 
a time interval with r S> At 2> thu The numerical data is 
compared to the analytical bound A| (solid line) given in 

Eq. (101. The fit (dashed line) suggests a (Dgjoo ~ 2~ 12K 
dependence. We average over at least 100 (40 for K = 10) 
sets of random coupling constants \4>k\ 2 and show the mean 
value with the standard deviation as blue points. Besides a 
magnetic held of b = 0.05 we use the totally mixed state for 
nuclear spins but polarize the electron spin in a direction. 


tity. We have also derived three different bounds for the 
time-averaged weak distinguishability and thereby recov¬ 
ered one previously known bound by means of a new 
method. 2 We have therefore been able to predict upper 
limits to the size of fluctuations in small closed quantum 
systems. 

Applying our analytical results to a QD setup in which 
an electron spin is coupled to nuclear spins of the host 
material through the hyperfine interaction enables us to 
make precise predictions. Since this spin system is typ¬ 
ically well isolated from its environment, QDs can be 
considered as a closed quantum system for sufficiently 
short time scales. We have simulated the time evolution 
of the total spin system and have analyzed its depen¬ 
dence on experimentally accessible parameters such as 
the strength of an external magnetic field and the polar¬ 
ization of the initial state of the nuclear spins. Intrigu- 
ingly, we have discovered cases in which strong magnetic 
fields do not prevent the electron spin from equilibration, 
while a polarized bath always diminishes the equilibra¬ 
tion capability. Finally, we have also investigated the im¬ 
portance of the number of nuclear spins on equilibration 
properties. We show both analytically and numerically 
that very small amounts of bath spins are sufficient to 
fully equilibrate the electron spin in our model. The an¬ 
alytical results even hold for a wider class of spin models, 
and, thus are not limited to our specific model. 


V. SUMMARY AND DISCUSSION 


In summary, we have shown how a general theory on 
equilibration can be applied to a realistic closed quantum 
system. We have introduced a specific understanding of 
equilibration relevant for our system under consideration 
and analyzed its properties by analytical calculations. 
Afterwards, we have applied this concept to a model of 
electron and nuclear spins in a solid state QD, which we 
have investigated by numerical simulations. 

A system is assumed to be in equilibrium, if an ob¬ 
server cannot distinguish, for most times, between the 
actual state of the system and its equilibrium state us¬ 
ing a finite set of measurements. Notably, two observers 
with different measurement sets could come to different 
conclusions. This equilibrium state is not necessarily a 
thermal state since it can, for instance, depend on the ini¬ 
tial state of the system. The distinguishability between 
the state of the system and its equilibrium state can be 
quantified by a suitable “measure of distance”. In this 
article, we consider the so-called weak distinguishability 
which vanishes whenever the system is in equilibrium. 
For two-outcome measurements, we have been able to 
show that a saturation of the corresponding expectation 
values is equivalent to equilibration. Furthermore, we 
have demonstrated how the variance of a time-dependent 
observable can be bounded by this weak distinguishabil¬ 
ity, which has allowed us to connect this abstract mathe¬ 
matical function to an experimentally measurable quan- 
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Appendix A: Saturation of variances 


In this Appendix, we show that one can separate 
(Dq) 4 = + A q in a time dependent part S^(t) that 

vanishes at large times and a time independent part Aq . 
To do so, we follow a previous analysis^ and use the fact 
that the matrix elements = {Ei\u>\Ej) of u> in energy 
space are given by 


UJij 


Pij F); — Ej 
0 else 


(Al) 








where Ei is the energy of the *-th eigenvector \Ei) of H. 
Now, we can rewrite the TAWD by 


(Do)t = j J (tr [p(t)6] - tr[AO]) 


= - / d i' 
t L 


(A2) 

2 


E^ 


-ih 1 (E„ — E m )t' \q 

nm'-' u j nm) K -'vnn 


t 


d t' 


'0 


_ w f e - i?L \E n -E m )t'Q 
/ J \h'nm '* L/ nm) z ' K - y n 


The last step is possible because for all n, m with 
w nm 0 follows E n = E m (see Eq. (All). There¬ 
fore, w nm e"^ £n_Em ^ is time independent. We define 
Va = (Pnm ^nm)Omn and gaps G a = E n E rn . Note 
that a is an abbreviation for a double index, running over 
all d 2 gaps. We then find 


1 


(D d )t = j / c W 


E' 


-i h- 1 G a t' 


(A3) 


= i f dt ' E^^ e ” m ” 1(G “” G ' 3K 

* 70 a,/3 

= iE( S + £ 

u a, 1 3 

G a =Gp G a ^C 

= ]Y, v ‘* v pJ* dt ' e ~ iH ~ 1(Ga ~ Gf<)t '+ E 


v a vte- in ~^- G ^' 


a,(3 

G a ^G p 


V a Vp . 

a,/3 
G a =Gp t^O 


=«aW 


—A. 


(A4) 


We can exclude the cases G Q = Gp = 0 in the sec¬ 
ond term because G a = 0 implies v a = 0. Note 
that S@(t) vanishes at least with j in limit of infinite 
times because the sum in 6 ( j{t) is upper-bounded by 
E{a,/3| G^Gf,} K v }\- Hence, 


lim (D d ) t = A a. 


(A5) 


Appendix B: Relation between distinguishability 
and variance 

Defining s(t) := \n(DQ) t , one can easily rewrite 
the definition of (D^) t (see Eq. Q) by 

D 6 (t) = [s(t) + l](D 6 ) t . (Bl) 

An average over the time interval I = [t,t + At] yields 

j dt' Dq( t') < ma x(Do)t' ■ f d t! [s(t') + 1]. (B2) 

Jl Jj 

With D d (t) = (tr[p(f)0] — tr[u)0]) 2 , the left hand side 
of the latter equation represents the variance of expec¬ 
tation values of O around the value tr[u)0] within the 


time interval I. Defining s(t) to be the average slope s(t) 
within J, we derive 

Var^i.Ai) < [m + l](D 6 ) u (B3) 

where we assume that ma x t i^i(D ( j)t' = (D^t- This as¬ 
sumption is correct if the system is on average approach¬ 
ing its equilibrium state. The value of s(t) is then nega¬ 
tive, however, s(t) > —1 holds strictly. This follows from 
both the semi-positive values of D^t) and the j in the 
definition of (D^t- Therefore, we prove that 

Var^A t)<(D d ) t . (B4) 

The latter bound holds for all systems that approach 
equilibrium in the sense defined above. If a system is 
already equilibrated, the TAWD (D^t is no longer de¬ 
creasing such that s(t) = 0. Note that in this limit, the 
estimate for the variance becomes exact. This is also the 
case at large times, where (Z?a) ( of each system saturates 
as we explain above in App. [A] 


Appendix C: Infinite time estimates 

In the following, we show how to estimate Ag using 
only basic information about the system. For this pur¬ 
pose, we start with the long-time limit of Eq. (A4) 


a 6 = E^ = EE ; 


7 7 

’>1 > 


(Cl) 


a,/3 

G a =Gf)yi 0 


j a,b 


where the sum in the last step is symmetrized by defin¬ 
ing a parameter j to run over all distinct values of en¬ 
ergy gaps, while a and b run over all rij gaps of size Gj. 
Therefore, v J a belongs to the a-th gap of size Gj ^ 0. We 
estimate the symmetric double sum by 


N 


N 


E^i = E W 2 + \ E ( XiX j + x i x i ) 

i i^j 

N 


N 


N N 

< e i^i 2 + 2 E (Ei 2 + N 2 ) 




N 


= ^E 


(C2) 


where {xi} is a set of N arbitrary complex numbers. Ap¬ 
plying this relation to Eq. (Cl), we obtain 


A<3 — E n i E 


(C3) 


which even is an equality as long as all gaps are not de¬ 
generate, i. e. rij = 1 V). With the maximum degeneracy 
of energy gaps Nq = max^ rij , we find 


A,<iVcEEKI 2 - 


(C4) 
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where both sums combined run over all d 2 gaps in the 
energy spectrum. In the previous notation this reads 

A 6 < Ar G ]TK| 2 - (C5) 


Third Estimate 


We derive estimate by following an approach along 
the lines of Ref. 0} We can then estimate Eq. (C6) by 


We now insert the definition of v a a nd use \p nm — w nm | < 
\pnm\i which follows from Eq. ( |A1| . Thus, we find 

Aq< N G Y,\PnmOmn\ 2 ■ (C6) 

n,m 


First and Second Estimate 

For the first estimate Ab, we use that \p n mO m n\ = 
\Pnm I I Omni and \p nm \ < ||/5||, where ||/5|| is the spectral 
norm of p. Using this and 0 = 0', this yields 

a 6 < A G ||p|| 2 ^o nm o; m 

nm 

= Nq ||/5|| 2 tr[0 2 ] 

=: Ab (C7) 

The same steps but p and 0 interchanged yields 

Aq < A^ := 1V G ||0|| 2 tr[p 2 ]. (C8) 


A q ^ Nq y ' PnnPmmOmnO71 


— A G trlpdiagOp^iagO], 


(C9) 


where \p n m\~ < PnnPmm because p is positive and 
(pdiag)nm = pnm&nm- With the Cauchy-Schwarz inequal¬ 
ity and tr[AR] < ||A|| tv[B] for positive A and B follows 


A 0 < Nq tr[0 2 p^ iag ] 

<iV G ||6|| 2 tr[p 2 iag ] 

<!V G ||6|| 2 tr[c5 2 ] 


N g 


IIOII 2 

deff 



(CIO) 

(Cll) 


This bound has previously been obtained in Ref. 0] on 
the basis of an analysis with pure initial states that have 
been expanded to mixed states afterwards. 
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